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ABSTRACT 

We discuss how the effective radius R e function (ERF) recently worked out by Bernardi 
et al. (2009) represents a new tcstbcd to improve the current understanding of Semi- 
analytic Models of Galaxy formation. In particular, we here show that a detailed 
hierarchical model of structure formation can broadly reproduce the correct peak in 
the size distribution of local early-type galaxies, although it significantly overpredicts 
the number of very compact and very large galaxies. This in turn is reflected in the 
predicted size-mass relation, much flatter than the observed one, due to too large 
(> 3 kpc) low-mass galaxies (< 10 11 Mq), and to a non-negligible fraction of compact 
(< 0.5 — 1 kpc) and massive galaxies (> 1O 11 M ). We also find that the latter discrep- 
ancy is smaller than previously claimed, and limited to only ultracompact (R e % 0-5 
kpc) galaxies when considering elliptical-dominated samples. We explore several causes 
behind these effects. We conclude that the former problem might be linked to the ini- 
tial conditions, given that large and low-mass galaxies are present at all epochs in 
the model. The survival of compact and massive galaxies might instead be linked to 
their very old ages and peculiar merger histories. Overall, knowledge of the galactic 
stellar mass and size distributions allows a better understanding of where and how to 
improve models. 
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1 INTRODUCTION 

The formation and evolution of early-type galaxies, 
characterized by a dominant central stellar bulge com- 
po nent, is still a ma tter of debate. The seminal paper 
by lEggen et al.l (|l962l ). postulated that stars are formed 
in a single burst of star formation from gas falling to- 
wards the center, and the evolution is passive thereafter. 
Such a simple scenario might be difficult to reconcile 
with the standard cosmological paradigm of structure 
formation, in which dark matter halos grow hierarchically 
through merging. The most advanced and up-to-date 
Semi- analytical mode l s (SAMs) of galaxy formation 
(e.g., ICole et all |2000|; iBenson et all 120031; iGranato et all 
2004 120061; iMenci et all 12004 ICirasuolo et all 120051 



Khochfar fc Burkertl [2003 : IVittorini et all 



2005 



Bower et al. 2006; Cattaneo et al.ll2006l: ICroton et all 


2006; 


Do Lucia et al. 


2006: Hopkins et al. 2006; Lapi et all 


2006; 


Shankar et all 


2006; Monaco et all 20071; ISomerville et all 



completely agree on the type of evolutio n undergone by 
massive galaxies (see also lDekel et all l2009l on the fraction 
of stellar mass formed in the initial, gas-rich burst of star 
formation, and on the role played by the late evolution 
driven by major and minor mergers. Nevertheless, all 
models agree that galaxies must have been much more 
compact at the epoch of formation, owing to a denser 
universe, larger gas fractions in the progenitors, and more 
dissipation. The latter prediction has been confirmed by a 
number of deep ob s ervations of high redshift galaxies (e.g.. 



Trujillo et al. 2006: 



van Dokkum et alj|200j; ICimatti et ail 



20081 ; ISaracco et all 20081 ). which have independently found 



early-type, high-redshift, massive galaxies to be a factor of 
a few more compact than local counterparts of the same 
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stellar mass. Note, however, that several observational bi- 
ases might limit the quality and reliability of some of these 
measurements fe.g.. iHopkins et al.l l2009dl : iMancini et al.l 
l2009l : Ivan Dokkum et aljhoogh . 



It is still debated how these compact galaxies have 
evolved from high-redshifts increasing their sizes in a way 
to fall on the size-mass relation we observe to d ay. A s 
already extensively discus sed by, e.g., IShen et all |2003), 
IShankar fc Bernardil (|2009h . Shankar et al. (2009a), the scat- 
ter around the local, median i? e -Af s tar relation decreases 
with increasing mass and, at high stellar masses, is nearly 
independent of the age of the galaxies, a challenge for most 
galaxy evolution models. In particular, older galaxies are 
observed to have a steeper size-mass relation than younger 
systems. 

One possible model put forward to explain the strong 
size evolution of the red, massive high-z galaxies is a 
sequence of minor, dry mergers. Such dynamical events 
can "puff-up" galaxies by adding mass in their out- 
skirts, efficiently increasing their sizes, although limiting 
the g rowth of the stellar mass within a factor of ~ 2 
/e.g.. iKhochfar fc Silk! l2006bl ; ICiottil 120091 ; iBernardil I200E 



- -jd3 — 

Bezanson et alJ 120091: [Cimattil I2OO9I: I Hopkins et alJ l2009~ 
Naab et alj 120091 ; Ivan der Wei et al.l 120091 '). Recent numer- 



ical simulations have however shed doubts on the actual 
efficiency of mergers in significantly puffing-up compact 
galaxies, and coherently bringing them along the rather 
tight structural re lations observed in the local Universe 
l|Nipoti et al l 120091 ). Moreover, hierarchical models suffer 
from the serious problem of failing in fully reproducin, 
the local size- luminosity relation fsee iGonzalez et al.ll200' 
Shankar et al. 2009a, and references therein). Another model 
proposed in the rece nt Literature to increase the sizes of 
massive galaxies is bv lFan et al.l l|2008h . They postulate that 
the evolution of compact galaxies undergoes a two-phase 
expansion: a first one caused by the sudden mass loss via 
quasar feedback, and a second one due to the slow mass loss 
via stellar winds during which the system slowly evolves to- 
wards a new equilibrium. Their model is broadly consistent 
with the data on the local size-mass relation. 

Given the large degree of freedom and significant uncer- 
tainties in current models of galaxy formation, it is necessary 
to look for other ways to test the validity of a given theory 
and/or discern ways to improve it. The aim of this paper 
is to provide such tests. We will show that the combined 
comparison with the size and mass distribution function of 
local galaxies can reveal interesting information on how to 
improve models of galaxy formation. In particular, in this 
work we will use a detailed hierarchical model of galaxy 
formation, show its failures and successes against available 
data, and discuss ways to improve it. 

We start in § [2] describing the data set we used. In § [3] 
we describe the hierarchical model adopted in this paper, 
and present its predicted size and mass distributions for 
spheroid-dominated galaxies. We will discuss in some detail 
the origin of the discrepancies between model predictions 
and the data, and possible ways to improve the model. We 
further discuss other aspects of the model in § 21 and con- 
clude in § [5] 



2 DATA 

The sample was extracted from the Data Release 4 of the 
Sloan Digital Sky Survey fSDSS; lYork et alj|200(t > with pa- 
rameters updated to the Data Release 6. The sample is mag- 
nitude limited to r-band deVaucouleur apparent magnitude 
14.5 < m r < 17.5. Spheroid-dominated objects were selected 
by requiring the concen tration index to be larger than 2.86 
|Nakamura et al]|2003l ). In the following we will take this 
dataset as our reference sample for early-type galaxies to 
compare the models with. This sample contains ~ 70, 000 
galaxies and it extends over a redshift range 0.013 < z < 
0.25, which corresponds to a maximum lookback time of 3 
Gyr. The effective radii and magnitudes computed by the 
SDSS reduction pipeline sy stem were corrected for sk y sub- 
traction problems following iHvde fc Bernard i (2009a). How- 
ever, as detailed in Bernardi et al. (2009), a sample selec- 
tion based solely on concentration, inevitably contains some 
SO and Sa galaxies. More caref ul cuts that provide clea ner 
samples are however possible l|Hvde fc Bernardil l2009al lbh . 
Therefore, in the following we w i ll also discuss the compar- 
ison with the IHvde fc Bernardil (|2009aT ) sample dominated 
by ellipticals. The comparison with the latter sample is par- 
ticularly meaningful for our purposes of understanding the 
evolution of the high-z, spheroid-dominated, compact, red 
and dead galaxies discussed in §[T] 



3 HIERARCHICAL MODELS OF GALAXY 
FORMATION 

3.1 OVERVIEW OF THE MODEL 

Bernardi et al. (2009) showed that the V/Vma,* method ap- 
plied to the samples discussed in § O yields the filled and 
open squares shown in Figure [4] for the SDSS early-type 
galaxies with concentr atioqj C r > 2.86 ( top pa nels) and 
those which satisfy the IHvde fc Bernardil (|2009al ) selection 
criteria (bottom panels) . The left panels of Figure [4] (dis- 
cussed in § I3.3[1 show the i? e functions (ERF hereafter), while 
the right panels the stellar mass functions (SMF hereafter) 
of the two samples. We defer the reader to Bernardi et al. 
(2009) for full details on the statistical and systematic un- 
certainties in computing a reliable estimate of the size and 
stellar mass functions of early-type galaxies. We just stress 
here that the V/Vmax method is one of the most appropri- 
ate and widely used techniques to compute the statistical 
distribution of galaxies of a given property (either luminos- 
ity, size, or vel ocity dispersion) a mong a flux-limited sample 
(see, e.g., also ISheth et al.ll2003h . For example, as detailed 
in Appendix[5] computing the ERF via some convolution of 
the luminosity function or the velocity dispersion function, 
provides consistent, although not as accurate, descriptions 
of the ERF. 

We compare our empiri cal determin a tion o f the ERF 
with the predictions of the iBower et al.l (120061 ) hierarchi- 
cal model (the 'Durham' model) . We will also b r iefly m en- 
tion some characteristics of the IDe Lucia et alJ (|2006l . the 



1 Defined to be the ratio of the scale which contains 90% of the 
Petrosian light in the r-band, to that which contains 50%. 



© 2009 RAS, MNRAS OOO.rflfHl 



3 



4 

3 
2 

1 

i.8 

0.8 
£ 0.6 
it 0.4 
0.2 

8:9 

0.6 




MINOR 



1 KLog M STAR <1 1.5 



DRY 
WET 



10.5<Log M STaR <1 1 JL 



10<Log M STAR <10.5 




MAJOR 



1 KLog M STAR <1 1.5 



10.5<Log M ST1R <1 1 J 




10<Log M STAR <10.5 



Figure 1. Comparison between the mean number of wet and dry mergers as a function of lookback time extracted from the merger tress 
of the Bower et al. (2006) catalog. Each row shows the merger history, averaged over 100 realizations, of galaxies with stellar mass at 
z = in three different mass bins, as labeled. In the left column we plot the mean number of minor mergers, with mass ratio < 1 : 3, while 
the right column shows the mean number of major mergers with mass ratio > 1 : 3. The dotted and solid lines refer to the mean number 
of wet and dry mergers, defined to have a (cold) gas-to-total mass fraction in the progenitors higher and lower than 0.15, respectively. 



'MPA model') model, although we will not compare di- 
rectly with their predictions as no available measure of the 
spheroid sizes are available from thei r catalogs at the mo- 
ment of writing. These models (see IParrv et alJ l2008l for 
a detailed description and comparison) follow the cosmo- 
logical co-evolution of dark matter halos, subhalos, galax- 
ies and supermassive black holes within the concordance 
ACDM cosmology. Both models have been implemented on 
top of the large, high-re solution cosmo logical N-body simu- 
lation MILLENNIUM RUN (|Springelll2005l ). The simulation fol- 
lows the evolution of N = 2160 3 dark matter particles 
of mass 8.6 x 10 8 ft -1 M , within a co-moving box of size 
500/i _1 Mpc on a side and force resolution of 5 h 1 kpc , 
from z — 127 to the present. The cosmological parameters 
adopted (fi m = 0.25, n b = 0.045, h = 0.73, fi A = 0.75, 
n = 1, and as = 0.9), are consistent with the combined 
analysis of th e 2dFGRS dColless et al.ll200lf ) and first year 
WM AP data (|Spergel et al.ll2003t ). The high mass and spa- 
tial resolution of the Millennium Simulation allows to track 
the motion of dark matter substructures inside massive ha- 
los, making it possible to construct merging history trees of 
all th e dark matter halos and subhalos inside the simulation 
box l|Sprin gcl 2005). The dynamical evolution of all satel- 
lite galaxies is followed until tidal truncation and stripping 
disrupt their host dark matter subhalos. Then a residual sur- 
vival time is estimated by computing the dynamical friction 
formula. 

To populate the dark matter subhalos with galaxies and 
black holes, both models adopt a set of equations to de- 
scribe the radiative cooling of gas, the star formation, metal 



enrichment and supernovae feedback, the growth and feed- 
back of supermassive black holes, the UV background reion- 
ization, and the effects of galaxy mergers. A full descrip- 
tion of the MPA and Durham models and a comparison be- 
t ween their main pred i ctions with ob s ervati o ns can be found 
inlBower et"afl (l200rJ).ICroton et all |2006ri.lDe Lucia et all 
(l2006t); |Pe Lucia fc BlaizotJ (|2007l ). iMalbon et all (|2007l ). 
iMarulli et al.l ( 20081 ). In the following, we will briefly sum- 
marize the main physical assumptions about the formation 
and evolution of spheroidal galaxies. 

In both models the morphology of a galaxy is deter- 
mined by the bulge-to-total ratio of its absolute, rest frame 
luminosity. Usually, a galaxy is classified as early type if it 
has Mbuigc — Mtotal < 0.4, where Mb u i gc and Mtotal are the 
B-band magnitude of the bulge and of the whole galaxy, 
respectively. The star formation is assumed to be propor- 
tional to the galaxy cold gas mass. Stellar winds and su- 
pernovae reheat a gas mass proportional to the mass of the 
stars, increase the metallicity of the disk material and in- 
ject substantial amounts of energy into their surroundings. 
Both models use the stel lar population synthesis model of 
iBruzual &: Charlotl (|2003l ') to describe the stellar population 
properties. 

The evolution of early type galaxies is then regulated 
by galaxy mergers, disk instabilities and AGN feedback. It 
is assumed that galaxy major mergers (Mgai^/Mga^i ^ 0.3) 
disrupt any stellar disk present and produce a spheroidal 
remnant, which contain all the stars of its progenitor galax- 
ies. The gas present in the merging galaxies forms stars in 
a burst, which are then added to the new spheroid. In the 
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case of a minor merger (M gsi i,2/M ga i,i < 0.3), the disk of 
the primary galaxy survives, and the stars and the gas of 
the satellite are added to the bulge and to the disk of the 
primary galaxy, respectively. While in the MPA model star- 
bursts occur in all minor mergers, in the Durham model 
there is a burst only if Mgaj^/Afgai.i > 0.1 and if the pri- 
mary galaxy has M gas /Mdi s k > 0.1. In both models AGN 
feedback in the radio mode further reduces or even stops 
late cooling flows in the halo centers. This behavior plays a 
key role in reproducing the exponential cut-off in the bright 
end of the local galaxy luminosity function, and in determin- 
ing the bulge-dominated morphologies and old stellar ages 
of the most massive galaxies in clusters. 

One key ingredient of hierarchical galaxy formation 
models of relevance for this paper, is how massive spheroids 
grow in size after they form. We here rec ap some of the 
basic model assumptions and refer to, e.g., .Gonzal ez et al.l 
(2009) for further details. The radius of the merger remnant 
is computed from energy conservation. Assuming virial equi- 
librium, the radius of the remnant is estimated solving the 
equation -E a fter = E\ + E2 + £"orb, with E\ and E2 the total 
internal energies of the two merging progenitors, and -Eorb 
their orbital energy. Each galaxy is characterized by its mass 
and radius, estimated assuming a R}^ n profile, with n = 4 
usually assumed for early-type galaxies. The energy of each 
galaxy is simply written as Et = kGMf / Ri, with Mi the 
total mass within its radius, k a constant depending on the 
profile, and Ri the half-mass radius. From the virial condi- 
tion, the radius of the remnant obeys the condition 



t Ml Mf . Ml M 2 ' V ' 

Ri + H 2 + z J° rb fli+iS 2 

with < /orb < 2 parameteri zing the (uncert ain) orbital 
energy of the progenitors (e.g.. lCole et alj|2000l ). 



3.2 THE GROWTH OF EARLY- TYPE 
GALAXIES IN THE MODEL 

One of the main drivers for the strong evolution of bulge- 
dominated galaxies in size and stellar mass in hierar- 
chical models is mergers (see also discussions in, e.g., 
De Lucia et~ai]|2006l ; |Pe Lucia fc BlaizotJ [20071 ; iParrv et all 
20081 ). Figure [T1 shows the mean number ./Vmergers of wet 
and dry mergers per Gyr a galaxy had since its formation 
epoch as a function of lookback time t. Each row shows the 
mean merger histor y, averaged over 10 realizations of the 
merger trees in the lBower et al.l ([2006:) catalog, of galaxies 
with stellar mass at z = residing in the mass bin indicated 
at the top of each panel, as labeled. In the left column we 
plot the mean number of minor mergers, with mass ratio 
< 1 : 3, while the right column shows the mean number 
of major mergers with mass ratio > 1 : 3. The dotted and 
solid lines refer to the mean number of wet and dry merg- 
ers defined to have a cold gas-to-total mass fraction in the 
progenitors higher and lower than 0.15, respectively. On av- 
erage, the number of dry mergers grows with final stellar 
mass, and galaxies that end up with M s t ar > 10 Mq tend 
to undergo significantly more numerous merging events than 
galaxies with lower final stellar mass. 

We can sketch a general trend for the evolution of mas- 
sive early-type galaxies in hierarchical models. A large por- 



tion of early-type galaxies is usually formed at high red- 
shifts through a wet merger of gas-rich disk progenitors and 
then continued accreting stellar mass mainly through minor, 
dry mergers. Therefore, the epoch of formation is generally 
identified to be at high-z, when the first starburst event 
took place and the central potential well was settled. We 
also find that a significant fraction of massive bulges is also 
formed through instability of gas-rich disk galaxies. Even in 
this case the initial size of the newly formed bulge is set 
by virial equilibrium and energy c onservation, adop ting a 
condition very similar to Eg. [T1 (see lCole et al.|[200oT for de- 
tails). Irrespective of their exact formation process, most of 
the massive spheroids in the model kept on accreting mass 
and increasing their sizes through dry mergers, larg ely ex- 
tending the epoch of stellar mass assembly. The lBower et al.l 
(2006) model predicts that all massive early-type galaxies, 
on average, undergo ~ 3 — 7 minor mergers and < 1 major 
dry mergers since their formation epoch (i.e., identified as 
the epoch of the major wet merger among the gas-rich disk 
progenitors). 

These late evolutionary features are a general trend for 
mos t hierarchical m o dels. F or example, we ha ve checked that 
the |Pe Lucia et all (|2006r i model (see also iGuo fc White! 
I2008T ) predicts a similar pattern for the growth of early- 
type galaxies, first characterized by a wet, formation phase 
peaked at high-redshifts, and a subsequent evolution dom- 
inated by a series of minor, dry mergers. The latter model 
similarly also predicts a rate of mergers increasing with in- 
creasing final stellar mass. 

The grey lines in the upper and lower panels of Fig- 
ure [5] show how the sizes and masses, respect ively, change 
for a h undred merger histories drawn from the lBower et al.l 
(2006) catalog, so that the final stellar mass at z = will 
fall within the mass bin labeled at the top of each panel. 
In each case the trees are followed back in time, choos- 
ing the most massive early-type progenitor, until this is no 
longer possible. We have checked that our randomly selected 
bulge-dominated galaxies in the model have mainly grown 
in mass through mergers, with disk instabilities increasing 
the spheroid masses by only < 10%, decreasing to < 2% 
for the most massive galaxies. This is clearly visible from 
the lower panels of Figure [2l where most of the stellar mass 
accretion histories remain almost flat most of the time and 
have sudden jumps in correspondence of merging events. 

Figure [2] shows that only a small fraction of low mass 
objects (with M sta r < 10 1 Mq) were early-types 10 Gyrs 
ago, but their mass has changed little since they first be- 
came early types (< 50%). The sizes of this population of 
galaxies, however, have grown by at least a factor of three 
over this time, as expected if their evolution is driven by mi- 
nor mergers. In contrast, a larger fraction of massive galax- 
ies had formed 10 Gyrs ago; however, whereas the mass of 
the population has grown by about a factor of three since 
then (with most of the growth occurring at lookback times 
greater than about 6 Gyrs) , the sizes have grown by a smaller 
factor (~ 2). Indeed, the models predict that a substantial 
fraction of massive galaxies (with M s t ar > lO 11 ]^©) have 
not changed in size for a long time. More specifically, in 
the highest mass bin considered, we find that about 21% of 
the galaxies have efficiently increased their sizes by a factor 
greater than 3, a significant fraction of 22% remains instead 
extremely compact (although still growing in mass), and the 
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Figure 2. Hundred merger histories for the evolution of the size (upper panels) and stellar mass (lower panel) as a function of time for 
galaxies of different mass at z = as predicted by the Bower et al. (2006) model. The median and variances of the merger histories are 
shown in each panel with red circles with error bars. 



rest undergoes a milder evolution with a size increase con- 
tained within a factor < 3. We have verified that while most 
(~ 60%) of the mass growth is added via major mergers, 
most (~ 90%) of the size growth, where this actually hap- 
pens, is via minor mergers. As extensively discussed in § 13.31 
the survival of a non- negligible fraction of compact and mas- 
sive galaxies until the present epoch is at variance with the 
SDSS distribution of early-type galaxies (with B/T > 0.5), 
and we will discuss possible causes and improvements to this 
problem in § 13.4.31 However, we will also show that when re- 
stricting the analysis to spheroid-dominated systems (with 
B/T> 0.7), the predicted number density of ellipticals with 
size ~ 1 kpc is consistent with the data. 



3.3 COMPARISON WITH THE SIZE AND 
STELLAR MASS DISTRIBUTIONS 

In this section we compare the data on the size and mass 
distributions of SDSS early-type galaxies with the predic- 
tions of SAMs. As partly explained above, it is beyond the 
scope of the current paper to build an ab initio complete hi- 
erarchical model able to predict the sizes, stellar masses and 
velocity dispersions of early-type galaxies, which we post- 
pone to future work. Instead, the aim here is to show the 



what we can learn from the simultaneous comparison with 
the size and mass distribution of early-type galaxies. 

3.3.1 THE PREDICTED SIZE-MA SS RELA TION 

The orange, green, and blue contour levels in Figure [3] in- 
dicate the region of the size-mass plane containing 68%, 
95%, and 99.7% of the whole SDSS sample of local early- 
type galaxies, respectively. All the contour levels have been 
wei ghted by the appropr i ate 1/ Knax, as in iBernardil (120091 ) 
and lShankar fc Bernardil l|2009l ). The cyan, yellow, and red 
contour levels contain the 99.7%, 95% and 68% of the cor- 
responding mock s ample of local early-t ype galaxies from 
iBower et all |2006). Figure [3] plots the iHvde fc Bernardil 
l|2009al ) sample, and therefore, for consistency, we only adopt 
mock galaxies with B/T > 0.7. We have also checked that 
the mean trend and scatter in the predicted size-mass rela- 
tion does not change significantly when adopting lower lim- 
its of the bulge component, such as B/T > 0.5. The solid 
squares with error bars show the predicted median and vari- 
ances of the sizes at fixed stellar mass for the models. We find 
that the model predicts an increasing size when moving from 
the lower to the highest masses up to M sta r~ 2 x 10 9 Mq. 
However, as seen in the Figure, at stellar masses above 
Af star ~ 2 x 10 9 M Q , at variance with the data the model 
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Figure 3. Predicted size-mass relation for the lBower et al.l j2006r ) 
sample, with the red, yellow, and blue contours marking the re- 
gion of the ii e -Mstar plane containing 68%, 95%, and 99.7% of the 
mock catalog, respectively, and the solid squares indicating the 
median R e at fixed stellar mass and its standard deviation. Also 
shown the observed size-mass relation for the SDSS sample with 
the orange, green, and blue contours marking the region of plane 
containing 68%, 95%, and 99.7% fraction of the total sample. The 
predicted size mass relation flattens at high-masses, at variance 
with the data, although still a significant comoving number den- 
sity of large galaxies is predicted by the model. The dot-dashed 
line in the figure is the best-fit relation from Shankar et al. 2009a, 
while the dotted lines roughly mark the locus of points 3-<r away 
from the median relation. 



predicts a strong flattening in the predicted size-mass rela- 
tion, while at M a t ar ~ 3 x lO n M0, the sizes start increasing 
again with increasing stellar mass. Similar fi ndings were also 
recently discussed bv lGonzalez et all <j2009h . 

Not only the observed and predicted size- mass distribu- 
tions differ in slope and zero point, but also the predicted 
scatter (in size at fixed stellar mass or vice versa), is much 
larger than the observed one. Figure [3] clearly highlights the 
problem in the predicted size- mass relation, alread y noticed 
by some previous studies (e.g., [Gonzalez ct al. 2009, Shankar 
et al. 2009a) . We will show below that despite some possible 
improvements towards reproducing the size and mass func- 
tions, the full match to the size-mass relation remains an 
extremely non-trivial task for the model. 



'.3.2 THE PREDICTED SIZE F UNCTION 



The blue area in Figure \4%i shows the iBower et all (|2006l ) 
predicted size distribution for early-type galaxies computed 
by counting the number of sources within a given bin of 
size R divided by the volume of the Millennium simulation, 
(500 h^ 1 Mpc on a side, where we set h = 0.7 to get the 
numb er density in units of Mpc 3 ). Following Gonzale z et al.l 
(2009), we further correct the half-mass radius R by a fac- 
tor of 1.35, assuming t hat all galaxies strict ly follow a de 
Vaucouleurs profile (see lGonzalez et al.ll2009l and references 
therein). We further divide the sizes by a factor of 0.7 to 
scale from the units of kpc in the model to kpc in the 
data. We also divide the masses in units of /i _1 Mq by 0.7 
and increased them by 0.05 dex to account for the slightly 
different initial mass functions adopted in the model (Ken- 



nicutt 1983) and the data (Chabrier 2003; se e Table 2 in 
Bernardi et al. 2009). We then select from the lBower et al.l 
(2006) catalog all those galaxies with a total stellar mass 
above logM sta r/AfQ ~ 9.5, consistent with our sample, and 
a prominent bulge component. The blue area in Figure [4]a 
includes all galaxies with Mb, bulge — Ms, tot < 0.4 (lower 
limit) and Mb, bulge — Mb, tot < 1.3 (upper limit), which cor- 
responds to the subsample of galaxies with a bu lge-to-total 
stellar mass ratio of BIT = 0.5 ± 0.2 (see, e.g., ICole et all 
|2000| . lLaurikainen et aL 20091 and refer e nces t herein) . As ex- 
tensively discussed in Gonzalez et alj (|2009l ). a cut in the 
mock catalog of B/T ^ 0.5 is equivalent to a cut in con- 
centration of C r > 2.86, the one adopted by Bernardi et al. 
(2009) to classify early-type galaxies in the sample shown 
with solid, filled squares in Figures b. 

Overall, the model predicts a nearly Gaussian-shaped 
size distribution, with a peak around ~ 3 kpc, which quite 
closely resembles the observed one. In particular, the model 
predicts a sufficiently large number of intermediate-size and 
large galaxies (Re> 1 kpc) in broad agreement with data. 
Such a result is not trivial. The outputs of hierarchical mod- 
els concerning the sizes of galaxies have been discussed sev- 
eral times in the Literature (see § [Q. However, emphasis 
has been usually put on the tendency of the models to pre- 
dict too many compact (R e < 1 kpc) galaxies, but not much 
was discussed about larger galaxies. We confirm indeed that, 
with respect to the data, the hierarchical model considered 
here (see solid line in Figure |4]i), shares this common fea- 
ture of predicting too many, up to a factor of ~ 10 higher, 
early-type galaxies with sizes below ~ 1 kpc and stellar mass 
M star > 3 x 1O 1O M0. On the other hand, it also predicts the 
number density of larger galaxies to be somewhat in bet- 
ter agreement with SDSS measurements, although it still 
overpredicts by a factor ~ 2 the number density of galax- 
ies with sizes between 7? e ~ 3 — 5 kpc, and even more that 
of ultra- large galaxies with R e > 15 kpc. These outcomes of 
the model would obviously not have been evident a priori by 
just plotting the size-mass relation, which predicts a rather 
flat distribution in the median sizes, and it lends support 
on the importance of comparing model outputs with size 
distributions as well as mass functions. 

Bernardi et al. (2009) showed that the CV > 2.86 sam- 
ple has a significant contribution from SO galaxies (> 30%). 
The latter class of galaxies has light profiles best-fitted by a 
combination of an exponential and a r -1 '' 4 de Vaucouleurs' 
profile (e.g.. Ide Vaucouleurs! IT959). which yields somewhat 
larger effective radii (by ~ 0.1 dex) than a pure r -1 / 4 pro- 
file. This in turn might bias the result of Figure |4]i, because 
by assuming a pure r -1 ^ 4 profile for all early-type galaxies, 
we could have overpredicted the number of compact galax- 
ies. We have however checked that even after increasing the 
sizes of all SO mock galaxies (those with B/T < 0.7) by 0.1 
dex to match the results of Bernardi et al. (2009), the pre- 
dicted ERF is very similar. We therefore conclude that the 
overproduction of early-type galaxies at small scales is real 
and inherent in the model itself. 

Interestingly, Figure shows that when comparing 
to the elliptical-dominated sample by iHvde fc Bernard! 
(2009aj), characterized by a minimal contribution from SO 
and Sa galaxies (open squares), the overproduction of com- 
pact galaxies is significantly reduced, and possibly confined 
to only ultracompact (R e < 0.5 kpc) galaxies. For this com- 
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Figure 4. Upper left panel: compar ison of the measure d ERF for galaxies with CV > 2.86 (filled squares) from Bernardi et al. (2009), 
compared to the one predicted by the Bower ct al. ( 2006) model. The different lines show the contribution of galaxies from the model with 
masses above a different threshold, as labeled, and with B/T > 0.5. The colored areas around each curve mark the uncertainty in defining 
an early-type galaxy in the model with the color regions bracketing the galaxies with B/T > 0.3 and B/T > 0.7. The red long-dashed 
line marks the predictions of a model in which the galaxies with sizes below and above the dotted lines in Figure [3] are removed from 
the mock catalog. Lower left panel: same format as panel a but now restring to the ERF computed from the Hyde & Bernardi (2009a) 
sample (open squares), dominated by ellipticals, compared to the predicted ERF for galaxies with B/T > 0.75. Right upper panel: 
comparison between the measured and predicted stellar mass functions for samples as in panel a. The data (filled squares) refer to the 
sample of galaxies with CV > 2.86. The solid line identifies the contribution of all early-type galaxies in the model with B/T > 0.5, while 
the red long-dashed line considers only B/T > 0.5 galaxies within the dotted lines in Figure [3] The other lines show the contribution of 
galaxies with effective r adius above a different thr eshold, as labeled. Right lower panel: same format as panel b but referred to the stellar 
mass function from the Hyde & Bernardi (2009a) sample (open squares), and for model galaxies with B/T > 0.75. The colored areas 
around each curve mark the uncertainty in defining an early-type galaxy in the model with the color regions bracketing the galaxies 
with B/T > 0.55 and B/T > 0.9. It is clear from these plots that simply removing the compact galaxies from the sample severely 
underpredicts the high-mass end of the stellar mass function. 



parison we used mock gal axies with B/T > 0.75, ro ughly 
the minimum B/T in the iHvde fe Bernardil |2009al ) sam- 
ple (solid line wit h blue area). Wh e n mov ing from the 
CV > 2.86 to the IHvde fc Bernardil |2009al) sample, the 
shapes of the observed and predicted ERF are similar: while 
the observed ERF only slightly decreases in normalization, 
the predicted size function significantly drops below < 1 kpc, 
better matching the data. 

iTruiillo et all (|2009l ) (see also iTavlor et alj 120091 ') have 
recently claimed that in the SDSS Data Release 6, only a 
tiny fraction (~ 0.03%) of compact (< 1.5 kpc) and massive 
(M s tar^ 8 x 10 10 Mq) sy s tems is present in the local Universe 
(z < 0.2). lTruiillo et al.l d2009) also claim that this estimate 



is actually an upper limit, given that most of these galaxies 
are metal-rich and relatively young (~ 2 Gyr). They there- 
fore conclude that the nearly complete absence of very com- 
pact and old galaxies in the nearby Universe, is at variance 
with predictions of detailed hierarchical models that grow 
galaxies through mergers. The res ult discussed in Figur e^, c 
are in broad agreement with the ITruiillo et al"] (|2009l ) find- 
ings, that hierarchical models, or at least the one considered 
here, tend to produce too many compact galaxies. However, 
we find that such a discrepancy is not as strong as claimed 
before, and we also add that the difference is morphology- 
dependent, i.e., significantly reducing with increasing B/T. 

On the other hand, when moving to more spheroid- 
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Figure 5. Same format as Figure [4] but here the model predictions have been modified by assuming that satellite galaxies with stellar 
mass below Mstar~ 3 X 1O 1O M0 have merged with the central galaxy. It is clear that increasing the rate of mergers of the less massive 
spheroids yields a considerably better agreement with the data, although it does not improve the match to the size-mass relation (see 
Figure [6} . 



dominated galaxies, the model tends to produce an higher 
fraction of galaxies with R e > 3 kpc, predicting up to a fac- 
tor of ~ 10 higher number density of super-large galaxies 
(with -Rc> 30 kpc). This overproduction of large galaxies 
with mass around ~ 10 M© is one of the major causes, 
together with the overproduction of compact and massive 
galaxies, for the flattening of the size-mass relation seen in 
Figure El 

3.3.3 THE PREDICTED STELLAR MASS FUNCTION 

Figures 0f>, d show the corresponding comparison between 
model predictions and SP SS data on the SMF for t he CV 
sample (panel b) and the iHvde fc Bernardil (|2009al ) sam- 
ple (panel d). It is apparent that the model provides a 
poorer match to the SMF, irrespective of the sample con- 
sidered. More noticeably, at variance with the data the 
model predicts a "bump" in the number density of early- 
type galaxies around M star ~ 1.5 x 10 10 M©, and falls short 
by a factor of ~ 2 in producing enough galaxies with 
ma ss M star > 10 11 Mq. We h ave checked that the model 
by |Pe Lucia fe Blaizoj ([2007) , possibly due to the differ- 
ent treatment of dynamical friction timescales (see, e.g., 



IParrv et all |200S| ; ISeek Kim etall l2009h . indeed produces 
a flatter distribution in the number density of early-type 
galaxies, in better agreement with the data. 



3.4 LOOKING INTO THE MODEL 

3.4.1 INCREA SING THE SA TEL LITE- CENTRA L 
MERGER RATE 

A very preliminary test to improve on the flatness of the 
predicted size-mass relation would be to simply remove from 
the mock catalog the large and less massive galaxies (R e > 3 
kpc and M a t ar ~ 1O 1O M0), and the compact and massive 
ones (-Re< 1 kpc and M star ~ 10 1 Mq). One simple way 
to do this would be to assume that those galaxies which 
are strong outliers with respect to the observed local size- 
mass relation in Figure [3] should be considered as irreg- 
ulars and/or disk-dominat ed galaxies instead of spheroids 
(e.g.. iHopkins et al.ll2009cl and references therein). The re- 
gion in the 7? e -M sta r plane within the dotted lines in Fig- 
ure [3] roughly defines the "region of acceptance" at a ~ 3a 
level. All the mock galaxies which lie below and above those 
lines can be safely considered as outliers. Recomputing the 
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Figure 6. Same format as Figure \3\ with S PSS data compared 
to the predicted size-mass relation for the iBower et all (2006) 
model with low-mass, type 1 galaxies merged with their centrals, 
as discussed in § 13,41 For comparison we also show with a long- 
dashed line the relation from Figure [3] The size-mass relation is 
still not properly reproduced. 



merger. More recently, several groups working with the same 
model considered here, have shown that an increased galaxy 
merger rate produces bette r fits to the luminosity fun ction 
for different Hub ble types ([Benson fc Devereuxl 120091 ) and 
for the clustering |Seek Kim et al. 20091 ). 

When merging the satellite galaxies, we have also up- 
dated the sizes of the centrals following Eq. {TJ and setting 
/orb = 2. However, as shown in Figure [5]ci, c, at variance 
with the SMF, the comparison to the observed ERF does 
not significantly improve. Note that using lower values of 
/ or b would yield larger sizes for the remnants, thus further 
worsening the match to the data. The model continues to 
predict too many large and ultra-large galaxies, with respect 
to those observed. More importantly, we find that the result- 
ing size-mass relation for this kind of model still presents 
similar discrepancies with respect to the data, as shown in 
Figure[6] Low-mass galaxies still have too large sizes, while a 
non- negligible fraction of the massive galaxies still have too 
compact sizes, thus preserving the flattening of the relation. 



predicted ERF and SMF only considering the subsample of 
galaxies within the dotted lines, yields the red long-dashed 
lines in Figures \^fi, b, c, d. It is clear that, irrespective of the 
actual physical basis motivating such drastic cut in the mock 
catalog of early-type galaxies, this is not a satisfactory solu- 
tion, as it severely further underpredicts the number density 
of massive galaxies, without significantly reducing the bump 
around M sta r~ 1.5 x 1O 1O M0, and also worsens the match 
to the ERF. 

We have checked that a significant fraction of the out- 
lier galaxies are satellites. Therefore, one possible way to 
remove the outlier galaxies might be to simply assume a 
higher rate of mergers with the central galaxy in the same 
parent halo. Figure [S]shows the model predictions modified 
by assuming that all satellite galaxies with stellar mass be- 
low M sta r= 3 x 1O 1O M0 have merged with their respective 
centrals. It is apparent how this hypothetical increase in the 
merger rate of the less massive spheroids yields a consider- 
ably better agreement with the observed SMF, lowering the 
number density at low masses, and increasing it at higher 
masses. 

Although the simple improvement proposed here for 
this model yields outputs in better agreement with the data, 
it requires a valid physical explanation and several indepen- 
dent tests (such as clustering). Some hints that may moti- 
vate a revision of the merger rates in the model have al- 
ready been discussed in the Literature. First of all, there 
is the observational fact that the frequently assumed minor 
merger rate of massive galaxies seems to be lower than what 
suggested by recent observations (Jogee et al. 2009). Sec- 
ond, although the galaxy merger rates broadly follow the 
high-resolution merger rates of the dark matter halos and 
subhalos in the Millennium simulation, the actual compu- 
tation of the former may still not be accurate enough. In 
fact, baryons can affect their s urrounding halos be cause of, 
e.g., adiabatic contraction (e.g.. iTissera et al 2009, and ref- 
erences therein), leading to a denser subhalo, which would 
not be as easily tidally stripped, hence allowing for a faster 



3.4.2 THE SIZE-MASS RELATION AT HIGHER 
REDSHIFTS 

Figure [7] shows the size- mass relation for ellipticals (B/T 
> 0.75) at different redshifts, as labeled. It is evident that 
the flattening discussed above is present at all epochs (see 
also Figure 5 in Shankar et al. (2009a)). At fixed stellar mass, 
galaxies with M ata r< 1O 11 M0, tend to shrink by a factor of 
~ 3 at higher redshifts. At higher masses, although a large 
fraction manages to grow (at z < 2) in size by a significant 
amount enough to saturate (and actually overproduce) the 
number counts of large galaxies in the local Universe, still 
a large portion of it remains compact, as we further discuss 
in the following section. 

We here note that, with respect to observations, low 
mass galaxies tend to be already quite larg e (R e > 1 kpc) 
close t o their formation epoch. As reviewed bv lMancini et al] 
(2009), most of the low-mass galaxies observed so far in deep 
surveys, have in fact on average a factor of < 2 smaller 
sizes than the ones predicted here at the same redshifts. 
This is due a combination of several processes. Some of the 
spheroids are generated by strong disk instabilities in the 
progenitor. The exact treatment of instabilities and of the 
energy conservation between the p rogenitor d i sk and the 
final bulge, h ave b een discussed by I Cole et ail l|2000r i and 
I Bower et al.l J2OO6) , and are still subject to significant un- 
certainties that might require some closer comparison with 
hydro-simulations to provide more reliable answers. The 
spheroids formed instead through the merging of disk domi- 
nated, gas-rich progenitors, could possibly be more compact 
than currently assumed in the present model, as more gas 
in the progenitors effectively produces more compact sizes 
than the ones estimated from dis sipationless virial rela tions 
of the type given in Eq. [T] (e.g., iHopkins et al.ll2009dl . and 
references therein). The second problem of the formation 
and survival of compact and massive galaxies is a problem 
of possibly different nature and we discuss it in some detail 
below. 
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Figure 7. Same format as Figure[3]where we plot the predicted size-mass relation of early-type galaxies (B/T > 0.7) for the lBower et al.l 
(2006) model at different redshifts, as labeled, and compared with SDSS data at z = for reference. The long-dashed line in each panel 
is the size- mass relation at z = 0, from Figure [3] The presence of a flattened relation at all epochs suggests that some of the problems 
might be linked to the initial conditions. 



3.4.3 COMPACT AND MASSIVE GALAXIES 

As discussed above, Figures^, 6 reveal that the model tends 
to predict a larger number of compact (< 1 kpc) and mas- 
sive (M sta r> (0.5 — 1) x W^Mq) early-type galaxies, with 
respect to the data. Such a result persists even after al- 
lowing for an increased rate of mergers for central-satellite 
galaxies required to improve the match to the observed SMF 
(Figure[5]3,&). The extreme compactness of massive galaxies 
plays a considerable role in the flattening of the predicted 
size-mass relation shown in Figure [3] We here attempt to 
explore possible reasons responsible for producing smaller 
sizes, or suppressing their later growth, in massive galaxies. 

First of all, we note that model galaxies tend to be 
significantly older than similar galaxies in SDSS. Figure [8] 
shows normalized age distributions for model (solid lines) 
and SDSS (long-dashed lines), massive and compact galax- 
ies, as labe led. Here the ages in SDSS are from the spectral 
analysis of lGallazzi et al. I (|2005l ), and, although they might 
be biased by a sys tematic 1.5 Gyr wit h respect to mass- 
weighted ages (e.g.. [jimenez et al.ll2007t ). they are still sys- 
tematically lower than the predicted ones. We recall that 
the analysis carried out bv lShankar fc Bernardil (|2009l ) and 
Shankar et al. (2009a), showed that, on statistical grounds, 
the trends of age with stellar mass and size do not signifi- 
cantly depend on the exact choice for the age estimator, al- 
though different methods might yield significantly different 
ages on an object-by-object basis. We find that while galax- 
ies in SDSS peak at ages around 5 — 8 Gyr (z ~ 1 — 1.5) and 
have a rather broad distribution with a significant number 
of older and younger galaxies (see the more detailed anal- 
ysis in Shankar et al. 2009a), model galaxies are, on aver- 
age, older, peaked around 11 — 12 Gyr (z ~ 3 — 5), and 
have a narrower distribution. This result holds irrespective 
of the exact bulge fraction cut considered in the model (i.e., 
the discrepancy is present for both SO and pure elliptical 
galaxies). The difference in formation epoch might be, at 
least in part, responsible for producing a larger number of 
extremely compact galaxies at low redshifts. The model in 
fact produces, at fixed stellar mass, more compact galaxies 
at higher redshifts (see Figure 5 in Shankar et al. 2009a), as 
naturally expected if older galaxies are born from gas-richer 
events. 

To further test the origin of the suppressed size evo- 



lution in a fraction of massive galaxies, we randomly se- 
lect three samples of compa ct, average, a nd la rge bulge- 
dominated galaxies from the iBower et all l|2006l) model at 
z — with stellar mass of ~ 1.5 x 10 11 M , and a me- 
dian size of R e =0A, 1.5, and 8 kpc, respectively. Figure [5] 
shows the comparison among the median Af a t ar and R e as a 
function of redshift (panels a and b, respectively) for three 
samples, showing that the compact galaxies start off as more 
massive galaxies at high redshifts but do not evolve, on aver- 
age, their original sizes, even if, surprisingly, they can grow 
in mass by up to a factor of ~ 10. In Figure^ we plot the 
overall median size-mass relation for the full set of progeni- 
tors of the galaxies examined in panels a and b. The galaxies 
included in Figure [9)3 are all the progenitors merging onto 
the main branch. In other words, we exclude both galax- 
ies along the main branch and those that merge with other 
galaxies that will subsequently merge onto the main branch. 
It is apparent that the compact galaxies have "progenitors" 
always characterized by sizes a factor of a few smaller than 
average, which we believe is one of the main causes for the 
negligible growth in their sizes. In fact, from the virial con- 
dition given in Eq. JTJ), it is clear that setting 77 = M2/M1 
and k = R2/R1, we get 

Rf = (1 + V) 2 /[l + V 2 /k + 2r?/ orb /(l + k)]R! . (2) 

Therefore, if the merging progenitors have comparable sizes, 
i.e., k ~ 1, then this translates into 7?/ ~ (1 + rf) 2 /[l + 
rf 2 + 2rf\ — Ri, irrespective of the value for 77, while pro- 
genitors with k » 1 will efficiently puff-up the size of the 
remnant by a factor ~ (1 + 77) 2 . Another possible concur- 
rent reason for the very compact remnants can be possibly 
associated with the choice of f or b in Eq. [1] Cosmological 
N-body simulations show in fact that merging dark matter 
halos mostly do so on app roximately parabolic orbits (e.g., 
iKhochfar fc BurkertfeOOq ). implying values of f or b ~ 0, and 
hence larger remnant sizes. 

Overall, we conclude that, wit hin the hierarchical 
framework of the IBower et al.l (|2006l ) model, allowing the 
compact, high-z galaxies to be born at more recent epochs 
and to possibly merge with more massive and larger pro- 
genitors, should in principle allow them to grow more in 
size yielding better agreement with the low portion of the 
observed ERF. We therefore confirm that the discrepancy 
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Figure 8. Normalized age distribution for the spheroids and bulges in the range of mass and size as labeled as predicted by the model 
(solid lines) and as calibrated in SDSS (long-dashed lines). There is a tendency for galaxies in the model to be older. This in turn might 
play a role in producing too many compact galaxies at low redshifts. 



between model predictions and data regarding the presence 
of compact and massive galaxies is real and it causes the 
flattening of the predicted size-mass relation. Nevertheless, 
the model can efficiently increase in size and mass a large 
fraction of spheroids, in a way to actually significantly over- 
produce the number of local, large galaxies, at least above 
> 10 kpc. 



4 DISCUSSION 

In the previous sections we showed that hierarchical models 
are characterized by a two-stage evolution, a fast wet, high- 
redshift phase, followed by a much longer assembly phase 
dominated by minor mergers (see Figure [TJ. A fraction 
(~ 10 — 50%) of the stellar mass of the spheroids is mainly 
formed through wet, high-z merging events among gas-rich, 
disk dominated subunits, while the rest of the stellar mass is 
added via later dry mergerfl This kind of evolution may be 
connected with the fast and slow accretion phases of single 
dark matte r halos seen in high resolution numerical simula- 
tions (e.g., Zhao et al. 2003; Lu ct al. 2006; Dicmand ct al. 
120071 : lAscasibar fc Gottloberll2008l ; IVass et al]|2009l 'l. These 
studies support a scenario in which during a first, fast, 
and chaotic phase, the halo builds up the central potential 
well, while during a second, much longer phase, matter is 
accreted in a smoother way. The latter phase could resemble 
a sequence of minor mergers, although the total amount 
of stellar mass carried towards the central regions during 
this l on ger phase is still un cl ear (e . g . , iBovlan-Kolchin et al.l 



20081; Kazantzidis et"aT1 120081; (Purcell et alj I200I 



Drorv fc Alvarez! 120081 ; Perez- Gonzalez et al.l 20081 . and 



references therein). An increasing stellar mass with time 
may be actually going on in disk-dominated galaxies, where 
the break radius in the light profiles has been observed to 
increase by a fact or of 1.3 since z ~ 1 (jAzzollini et aljfeoosl ; 
iBakos et al1l200Sl ). 



2 We stress that the growth of the bulge stellar mass via disk 
instability in the massive spheroids considered here is not more 
than 10%, and only ~ 2% for the galaxies with final stellar mass 
Af st ar> 1O U M . 



Minor mergers might also be at t he origin of the 
boxy and disky early-type dichotomy (e.g.,|C appcllar i et al 



2007 ; lEmscllc m et al ] |2007l ; lKang' et a . 2007; |Pasquau' et a 



20071 . and references therein). Assuming that wet and 



dry mergers are actually responsible for determining the 
disky and boxy nature of stellar orbits in ellipticals, 
then the results in Figure [T] would imply that older 
galaxies, which have undergone more dry mergers since 
their formation epoch, would end up, on average, with 
boxier isophotes with respect to younger galax ies (e.g., 
iKhochfar fc Burkertj[2006l; IKhochfar fe Silk|[2006al '). The re- 
cent study by Kormendv et al.l (120081) actual l y confi rms that 
older galaxies appear boxier. Almeida et al.l 1 20071 ') also dis- 
cussed that several other correlations, including the veloc- 
ity dispersion, size and luminosity are reproduced by the 
I Bower et al.l (|2006t ) model with no extra tuning of the pa- 
rameters. 

During the wet phase, the correlations between 
the central black hole mass Mbh and their host 
galaxy potential wells, characterized by their a, might 
have also been settled, especially in A GN feedback- 
constrained galaxy evolution models (e.g., Granato et al.l 
|2004 iHopkins et al.ll2006l; iMonaco et al.ll2007l ; but see also, 
e.g.. lMiralda-Escude fc Kollm eicr 2005). If most of the black 
hole mass was already in place at the end of the wet phase, 
but only about half of the host galaxy mass was assembled, 
the hierarchical models considered here would then natu- 
rally imply an higher black hole mass to stellar mass ratio 
with respect to the local one, i.e., a positive evolution in 
the normalization of the MBH-Af ata r relation. On the other 
hand, given that the host halo potential well is rapidly built 
during the fast accretion phase, it is also reasonable to ex- 
pect the Mbh-o" relation to possibly already be fully estab- 
lished at the epoch of the wet phase (e.g., Granato et al.l 
12004 iMarulli et al.ll2008l ; IHopkins et atll2009bl ). Some em- 
pirical works have, in fact, found only marginal evidence 
for evolution in the Mbh-o" rel ation, and possibly only in 
the more massive systems, ( e.g., IShields et al.ll2006l ; lGaskelll 
2009; Sha nkar et al.l l2009bl . and references therein). How- 
ever, minor dry mergers are expected to ha ve some im- 
pact on the initi al velocity dispersion a (e.g.. Iciottil 120091 ; 
iNaab et al]|2009f ). Also, late black hole re-activations (e.g., 
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Figure 9. Median M B t ar and R e as a function of rcdshift (panels a and b, respectively) for three samples of galaxies randomly extracted 
from the I Bower et al . (2006) catalog with similar stellar mass at z = but very different final sizes (compact, large, and average, as 
labeled). Panel c shows the median radius versus stellar mass for the progenitors of the three samples. The progenitors of the compact, 
massive galaxies systematically have lower sizes. 



Menci et al. 2004, Vittorini et al. 2005), and/or black hole 
mergers (e.g., Volonteri et al. 2005, Malbon et al. 2007), 
might have increased the black hole masses since the wet 
epoch, further influencing the evolution in the scaling rela- 
tions between black holes and their host galaxies. Overall, 
the dynamical evolution of galaxies and their central black 
holes, tested against the local velocity dispersion function 
(Sheth et al. 2003, Bernardi et al. 2009) and fundamental 
plane of early-type galaxies, should provide valuable addi- 
tional insights into our understanding of galaxy evolution. 



5 CONCLUSIONS 

In this paper we make use of the data sets derived from 
SDSS DR6 by Bernardi et al. (2009) and lHvde fc Bernardil 
l|2009al ), used to derive the size and stellar mass functions 
for a sample of early-type galaxies with concentration C r > 
2.86, comprised of both ellipticals and SO galaxies, and a 
sample dominated by ellipticals, respectively. 

We compare the se statistical distrib utions with the hi- 
erarchical model by iBower et al l (|2006t ). The aim of this 
exercise is to show how the simultaneous comparison of the 
size and mass distributions can reveal interesting insights 
on how to improve the performance of theoretical models of 
galaxy evolution. 

We find, in agreement with previous studies, that this 



hierarchical model provides a poor match to the size-mass 
relation of local galaxies, irrespective of the exact sample 
we compare it with. In particular, the model tends to pro- 
duce a much flatter relation than the one actually observed. 
This flattening is mainly produced by the combined effects 
of having, with respect to the local data, too large (~ 3 
kpc) low-mass galaxies (< 1O 11 M0), and of having a non- 
negligible fraction of compact galaxies (< 0.5 — 1 kpc) at 
high masses (> lO n M ). 

Such discrepancies are reflected in the predicted size dis- 
tribution. Although the model produces a size distribution 
in broad agreement with the data, it tends to overproduce 
the number of large galaxies beyond the peak (> 3 kpc), 
and the number of very compact galaxies (< 1 kpc). We 
discussed that the former issue is present at all epochs, and 
it might therefore be linked to how spheroids are formed 
in the first place, either from not properly treating initial 
disk instabilities and /or computing the sizes of remnants in 
gas-rich mergers. 

Regarding the overproduction of compact and massive 
(M st ar~ (0.5 — 1) x 10 11 M Q ) galaxies with respect to the 
data, already pointed out in the recent Literature, we find it 
to be less prominent than previously claimed, and confined 
to only ultracompact galaxies (J? e < 0.5 kpc) when consid- 
ering only ellipticals. We discuss two possible reasons be- 
hind the survival of such compact galaxies until the present 
epoch. First, we find that model early- type galaxies tend to 
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be significantly older than those in SDSS. This in turn might 
induce more compact galaxies at fixed stellar mass, given 
that galaxies formed at higher redshifts are more compact 
(see Shankar et al. 2009a). We also find that model early- 
type compact galaxies underwent peculiar merging histories 
characterized by extremely compact progenitors, that could 
prevent them to efficiently grow their sizes. 
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APPENDIX A: PREDICTING THE SIZE 
FUNCTION FROM CONVOLUTION METHODS 

The filled squares in Figure Al represent the Bernardi et 
al. (2009) estimate of the ERF obtained from the V/V max 
method for our sample of early-type galaxies selected with 
concentration C r > 2.86. 

For consistency, we here show that the VVKnax-based 
ERF is consistent, within the errors, with the one obtained 
from the convolution of the luminosity function or velocity 
dispersion function with the bivariate distribution of points 
in the L — R e plane. Mor e speci fically, following the method s 
outlined in ISheth et all (120031 ) and IShankar etall l|2004t) . 
we have convolved the luminosity function <&(!/) with the 
bivariate distribution of Lj and R e , 

HK) (Al) 

3 

Here £ij is the fraction of sources in the sample with effec- 
tive radius R\ and luminosity Lj, normalized to the total 
number of sources with luminosity Lj . The luminosity func- 
tion $(1/) has been computed from the VVVmax method by 
Bernardi et al. (2009) for the same sample of galaxies, and 
we refer the reader to that paper for analytical fits and de- 
tailed discussions of the sample. The result of Eq. (|A1|) is 
shown in Figure Al with long-dashed lines, which bracket 
the statistical uncertainties in the luminosity function fit 
parameters. We have also used the bivariate distribution in 
equation (|A1|) applied to the velocity dispersion function 
$(cr), again derived by Bernardi et al. (2009) for this same 
sample, and with the weights £ij now computed from the 
distribution of sources in the a — R e plane. The result is 
shown with dotted lines in the same Figure, again brack- 
eting the statistical uncertainties in the velocity dispersion 
function fit parameters. 

This exercise proves that, as expected, convolutions of 
other statistical distributions $(:r) with their appropriate 
scaling relations x — R e , provide consistent results. How- 
ever, it also shows that the accuracy of the results relies on 
the accuracy of the input weights &(x), that in turn proves 
the importance of directly adopting the V/V max method to 
derive a more precise estimate of the ERF. 
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Figure Al. Size function of the CV > 2.86 sample obtained from the V/Vmax method (solid squares), compared with the size functions 
obtained from the convolution of the luminosity function (long-dashed lines) and velocity dispersion function (dotted lines), with the 
bivariate distributions of points in the L — R e and a — R e planes, respectively (see text). 
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